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Abstract 

Equilibrium eigenstates of an axisymmetric magnetically confined plasma with toroidal 
flow are investigated by means of exact solutions of the ideal magnetohydrodynamic 
equations. The study includes "compressible" flows with constant temperature, T(^), 
but varying density on magnetic surfaces and incompressible ones with constant density, 
g(ip), but varying temperature thereon (the function ip(R,z) is the poloidal magnetic 
flux-function with (R,z,4>) cylindrical coordinates). These variations are necessary for 
the existence of tokamak steady states with flow. The "compressible" and incompressible 
solutions are associated with the ansatz u 2 /T oc M = constant and d(gu> 2 )/dip = Ai/j, 
respectively, where ui{ip) is the rotation frequency. In both cases eigenfunctions of the 
form ip£ n = Zi(z)TZ n (R) (£, n = 1, 2, . . .) describe configurations with £xn magnetic axes. 
Owing to the flow, the respective eigenvalues can be considered in two alternative ways: 
either as flow eigenvalues, Mq and A e , which depend on a pressure parameter, Pq, or as 
pressure eigenvalues, Pq, which depend on the flow parameters Mo or A. In the context 
of the latter consideration when the flow parameters are varied continuously there are 
transition points, (Mo) m and A m (m = 1, 2, . . .), at which an additional magnetic axis is 
formed. This flow-caused change in magnetic topology is possible solely in the presence 
of toroidicity because in the limit of infinite aspect ratio the axial flow does not appear in 
the equilibrium equation. Also, the lower the aspect ratio the smaller (M ) m and A m . In 
addition, the effects of the flow and the aspect ratio on the Shafranov shift are evaluated. 



I. Introduction 

Over the last decades it has been established experimentally and theoretically that the 
flow affects the confinement properties of magnetically confined plasmas. In particular, 
the flow and especially the flow shear play a role in the formation of edge transport 
barriers (L-H transition) as well as of Internal Transport Barriers (ITBs), two enhanced 
confinement modes in tokamaks (e.g. see Ref. pQ and Refs. cited therein). Also, the 
majority of the advanced tokamak scenarios include flow. The ITBs usually are associated 
with reversed magnetic shear profiles. In addition, a possible magnetic topology of static 
(no-flow) equilibria with reversed current density in the core region proposed in Ref. [2] 
consists of multitoroidal configurations having non-nested magnetic surfaces. 
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Magnetohydrodynamic (MHD) equilibrium equations for axisymmetric magnetically 
confined plasmas with (a) isothermal magnetic surfaces and toroidal flow and (b) incom- 
pressible flow of arbitrary direction have been obtained in Ref. |3] and jlj, respectively. 
Although purely toroidal axisymmetric flows are inherent incompressible because of sym- 
metry, the former equilibria can be regarded as "compressible" in the sense that, alike 
the latter ones, the density varies on magnetic surfaces. Respective exact solutions were 
constructed in Refs. [HI E| and jH [Tj |H] and the impact of flow on certain equilibrium 
characteristics was examined therein. In particular, in. Ref. [7j we extended the well 
known Solovev solution EH] to unbounded incompressible plasmas and found that the 
flow and its shear can change the magnetic topology thus resulting in a variety of novel 
configurations of astrophysical and laboratory concern. 

The aim of the present study is to examine the possible impact of the flow on the 
magnetic topology for equilibria relevant to plasmas of fusion devices. The study includes 
equilibria with both "compressible" and incompressible toroidal flows. The main conclu- 
sion is that the flow in conjunction with toroidicity can change the magnetic structure 
by the formation of addition magnetic axes. The role of the toroidicity is important, i.e. 
this formation is not possible in the limit of infinite aspect ratio. 

The outline of the report is as follows. A derivation of the " compressible" and incom- 
pressible equilibrium equations in a unified manner is first reviewed and respective exact 
solutions are presented in Sec. II. In Sec. Ill equilibrium eigenstates of a magnetically 
confined plasma surrounded by a boundary of rectangular cross-section and arbitrary as- 
pect ratio are constructed. On the basis of these eigenstates the impact of the flow on the 
magnetic topology is then studied in conjunction with the role of toroidicity. The effects 
of the flow and the aspect ratio on the Shafranov shift are examined in Sec. IV. Finally 
the conclusions are summarized in section V. 



II. Equilibrium equations and solutions 



The ideal axisymmetric MHD equilibrium equations for the cases of "compressible" 
and incompressible toroidal flows are reviewed in this Section. In particular, a unified 
derivation is given without adopting from the beginning relevant energy equations or 
equations of states. They will be specified when necessary later. This rather detailed pre- 
sentation aims at making the discussion in the subsequent sections tangible, particularly 
as concerns the role of toroidicity. 

The starting equations written in standard notation and convenient units are the 
following: 

V ■ (gv) = (1) 
p(vV)v = JxB- VP (2) 

V x E = (3) 

V x B = J (4) 
V-B = (5) 

E + vxB = (6) 
An energy equation or equation of state (7) 

For axisymmetric magnetically confined plasmas with toroidal flow the divergence-free 
magnetic field and mass flow can be written in terms of the scalar functions ip(R,z), 
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I(R, z) and K(R, z) as 



B = JV0 + V0 x 

£ V = KV(p. 

The toroidal current density is then put, by Ampere's law, in the form 

J = A*^V0 - V0 x V/ 



(9) 



(10) 



Here (R, <p, z) are cylindrical coordinates with z corresponding to the axis of symmetry, 
ip labels the magnetic surfaces and A* is the elliptic operator defined as P 2 V • (V/P 2 ). 

By projecting the momentum equation (J2J) and Ohm's law © along the toroidal 
direction, the magnetic field, and perpendicular to the magnetic surfaces some integrals 
are identified as flux functions, i.e. functions constant on magnetic surfaces, and Eqs. 
(JU-6) are reduced to simpler ones. In particular, the V0-component of (j2J) yields 



V0 • (W x V/) = 



111 



implying that I = I(if>). Also, expressing the electric field in terms of the electrostatic 
potential, E = — V$, the component of Ohm's law along B leads to 



B ■ V$ = 0. 



(12) 



Eq. ()12j) implies that $ = <&(ip), viz. E is perpendicular to magnetic surfaces. One 
additional flux function is obtained by the projection of Ohm's law along Vip: 



'dip qR 2 ^ ^ ^ 



Therefore, the quantity 



K 



0J 



dip ' 



(13) 



(14) 



identified as the rotation frequency, is a flux function uo = uj(ip). 

With the aid of equations ((TTJ)-([T3J) the components of equation (J2J along B and 
respectively yield 



VP 



u 2 R 2 



[A*iP + II']\ViIj\ 2 + R 2 



VP - £cj 2 v( 



2 

P 2 

T 



•B = 

V^ = 0, 



(15) 
(16) 



where the prime denotes differentiation with respect to ip. 

In order to reduce Eqs. (JT5j) and (|TB|) further an energy equation or an equation of 
state is necessary. Owing to the large heat conduction along B, isothermal magnetic 
surfaces, T = T(ip), is an appropriate equation of state for fusion plasmas. In this case 
employing the ideal gas law, P = XgT, integration of (|T5j) yields 



P = P^)exp( 



u 2 R 2 
2XT 



(17) 
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where P s (ip) is the pressure in the absence of flow. In the presence of flow the pressure 
and therefore for T = T(ip) the density are in general not constant on magnetic surfaces, 
thus giving rise to " compressibility" . 

With the aid of (JT7J), Eq. (JTHj) leads to the final "compressible" equation 



A* if) + II' + R 2 



' ' 2 \XTJ 



exp( w )=0. 



Eq. (fTH|) was originally obtained in Ref. 

An alternative equation of state is incompressibility: 

V-v = 0. (19) 

Consequently, (JTJ) implies that the density is a flux function, g = g(ip), and integration of 
(ITK1) yields 

P = PsW + ^- (20) 



Eq. (JTHJ) then reduces to 



A*?p + II' + R 2 P' S + ^(gtu 2 )' = 0. (21) 



This is a particular form of the axisymmetric equilibrium equation for incompressible flow 
of arbitrary direction obtained in Ref. jl]. 

Both equations ()18|) and (|21|) contain four flux-functions, three out of which, i.e. P s , 
I and uj, are common. The fourth function is T for the "compressible" equation and g for 
the incompressible one. For vanishing flow (fTH|) and (|2*T|) reduce to the Grad-Shafranov 
equation. The flow term in (fTH|l depends on uj and g through uj 2 / XT and its ^-derivative 
(shear) while solely the shear of the flow term guj 2 appears in (|2ip. 

Linearized forms of Eqs. (JTSj) and (|21|) in connection with appropriate assignments of 
the free flux functions they contain can be solved analytically. In the present study we 
will employ exact solutions as follows. 

"Compressible" flow 

The ansatz used to linearize Eq. (fTHj) is [3] jB] 

I 2 = Iq + Ify 2 

P s = 2P ^ 2 (22) 

u 2 1 M 2 



XT R 2 



constant 



Here, Io/Ris the toroidal vacuum field, the parameter I\ describes the magnetic properties 
of the plasma; P , 7, and M are a pressure parameter, the ratio of specific heats, and 
the Mach number with respect to the sound speed at a reference point (z — 0, R — Ro) 
with R to be specified later. Note that the toroidal current density profile can vanish on 
the plasma boundary via (jHJ). 

Eq. (fTHj) then has a separable solution, H(R)Z(z), when the constant of separation is 
equal to I±. For configurations symmetric with respect to mid-plane z = this solution 
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is written in the form 



ip(x,y) = Ci 



Jo 



2t\ e 



7 M 



+ C 2 Y 



7 Mn^ x'2 



cos (hy), 



(23) 



where x = R/Ro and y = z/R ; J and Y are zeroth-order Bessel functions of first- and 
second-kind, respectively; and r 2 = AP Q Ri. 



Incompressible flow 

In this case the ansatz employed to linearize (J21j) is jS] 

P = J 2 + J^ 2 

P s = 2P ^ 2 
K 2 v 



(24) 



gR 4 



2Aip 



A remarkable difference of this choice compared to (}22|) is that the flow term has non-zero 
shear \{qoj 2 )' ^ 0] and consequently A can be positive or negative. Also, note that, unlike 
Mq in ()22|). A is dimensional. 

A separable solution is now expressed in terms of the first- and second-kind Airy 
functions, At and B iy as |Hj 



i/)(x,y) = Ci 
+C 2 5i 



ARo\- 2 / 3 /AR: 
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(^ 3 -p,flJ)) 

cos (hy). 



(25) 



III. Multitoroidal eigenstates associated with the flow 



In connection with solutions (}23|) and (}25|) we are interested in the steady states of a 
tokamak the plasma of which is bounded by a conducting wall of rectangular cross-section, 
as shown in Fig. ^ In addition, we assume that the plasma boundary coincides with 
the outermost magnetic surface; accordingly, the function ip should satisfy the following 
boundary conditions 

Hy±) = (26) 

and 

</>(x ± ) = 0, (27) 

where y± = ±a/i?o an d x± = 1 ± b/Ro. The equilibrium becomes then a boundary- value 
problem. Eigenstates can be determined by imposing conditions ()26|) and (|27p directly to 
solutions (J23|l and (|23j) . Specifically, (|2Tj) applied to the z-dependent part of the solutions 
yields the eigenvalues 

J i = ~(^-f), «=1,2,... (28) 

for the quantity 1% which is related to the poloidal current function The respective 

eigenfunctions are associated with configurations possessing £ magnetic axes parallel to 
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Figure 1: The cross section of the plasma boundary. The aspect ratio a is defined as 
Ro/b. 



the axis of symmetry. Condition ([27]) is pertinent to the independent part of the solution. 
Owing to the flow this part contains the parameter Mq in the compressible case and A 
in the incompressible one in addition to the pressure parameter Pq. (To make further 
discussion easier we introduce the symbol F representing either M or A; this is particu- 
larly convenient to formulate results which are independent of "compressibility".) Thus, 
condition l2*7j) can determine eigenvalues depending on the parameter Pq which remains 
free, F n (Po) (n = 1, 2, 3, . . .), or vice versa, pressure eigenvalues Pq{F) with F being free. 
This parametric dependence makes the spectrum of the eigenvalues broader in comparison 
with the static one. The eigenvalues depend also on the geometrical quantities Rq and b 
but not on a (see Fig. [TJ); thus, the results to follow are independent of elongation a/6. 
The other parameters C\ and C<i contained in (|21j) and (J23j) are adapted to normalize ip 
with respect to the magnetic axis and to satisfy the boundary condition (}2Tj) respectively. 
Also, in the rest of the report dimensionless values of A will be given normalized with 
respect to 1 Kgr/(m 7 T 2 s 2 ). The eigenvalues F n (P ) and Pq(F) can be calculated nu- 
merically. For a fixed value of Pq, F n satisfy for all n the inequality F n+1 > F n . A similar 
relation is satisfied by Pq for a given value of F. The respective eigenf unctions are con- 
nected to configurations having n magnetic axes perpendicular to the axis of symmetry. 
Therefore, the total equilibrium eigenf unctions ipe n = Zi{z)lZ n {R) describe multitoroidal 
configurations having £ x n magnetic axes. For example, a static doublet configuration, 
-012, possessing two magnetic axes parallel to z-axis was studied in Ref. [TTj . 

Henceforth, for the sake of simplicity we will restrict the study to eigenfunctions ipi n 
describing multitoroidal configurations with n magnetic axes along the mid-plane z = 0. 
As an example the ^12-configuration for "compressible" flow is shown in (Fig. |2J). For 
vanishing flow, tokamak multitoroidal configurations of this kind were investigated in Ref. 
[2] in connection with hole current density profiles which can reverse in the core region; 
current reversal is possible because the configurations have non-nested magnetic surfaces. 
It is noted that for static equilibria only pressure eigenvalues are possible. In the presence 
of flow we examined eigenstates with the lowest of the pressure eigenvalues Pq(F) by 
varying continuously the flow parameter F starting from a value close to the first-order 
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Figure 2: The double-magnetic- axis toroidal "compressible" configuration in connection 
with the eigenfunction ip 12 with eigenvalue Mq = 1.692 for aspect ratio a = 3. 



static one, Pq{F « 0). It should be noted here that for incompressible flow, A = does 
not necessarily imply static equilibrium because of the non-zero shear [see Eq. (J24J) ]. It 
turns out that there are transition points F m (m = 1,2, . . .) at which the configuration 
changes topology by the formation of an additional magnetic axis (The subscript m here 
indicating a transition point must not be mixed up with the superscript n indicating the 
order of an eigenvalue). For incompressible flow this is shown in Fig. El Specifically, the 
singly toroidal configuration |3fa) has eigenfunction ipu with eigenvalue Pq(A = 0.09). 
By varying the flow parameter the pressure (eigenvalue) decreases and the configuration 
is shifted outwards and is compressed in the outer region while the pressure becomes 
lower (jSl^b)). Then, as the flow parameter reaches the first transition point Ax = —0.01 
a second magnetic axis is formed in the outer region close to the boundary causing the 
first one to be shifted inwards (JSfc)). This transition point is associated with pressure 
eigenvalue Pq{A = A\). By a further increase of the flow (lower negative values of A) the 
outer magnetic island becomes bigger and the whole configuration is shifted outwards up 
to the next transition point A 2 = —0.1 associated with pressure eigenvalue P$ (A = A 2 ) 
(Fig. Efd)) at which as before a third magnetic axis is formed in the outer region; thus, 
the configuration consists of three magnetic islands (Fig. Efe)) and is shifted outwards as 
the flow increases even more (Fig. E[f)). This procedure is carried on until the formation 
of n magnetic axes and it is also possible for "compressible" flow. It can approximately 
be regarded as a quasi-static "evolution" of the plasma due to the flow variation through 
flow-depended pressure equilibrium eigenstates. Therefore, it is not possible for static 
equilibria. An animation showing this "evolution" is available in Ref [T2]. Alternatively, 
varying continuously the pressure parameter P one can find pressure transition points, 
(Po)m, associated with flow eigenvalues. 
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Figure 3: A sequence of graphs showing the quasi-static "evolution" of the configuration 
as the incompressible flow parameter A is decreased for aspect ratio a — 3. The values of 
A for the individual plots are (a) A = 0.09, (b) A = -0.01, (c) A = -0.02, (d) A = -0.1, 
(e) A = -0.122, (f) A = -0.2. 



For aspect ratio a = 3 and compressible flow the first transition velocity from a singlet 
to a doublet configuration is of the order of 10 5 m/s. Velocities of this order have been 
measured in tokamaks and therefore the change in magnetic topology by the flow can be 
experimentally examined. 

It is emphasized that the change of the magnetic topology due to the flow is not 
possible in the limit of infinite aspect ratio because in this limit the equilibrium equations 
do not contain the axial velocity regardless of " compressibility" . Indeed, for a cylindrical 
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plasma of arbitrary cross-section the equations respective to (JT5j) and ((TB)) read 

B • VP = (29) 

V 2 ^ + ( p + : y)' = °- ( 3 °) 

Eqs. (|29|) and (J30j) follow respectively from (16) and (17) of Ref |T3] for vanishing poloidal 
velocity (F' = 0). Also, note that the pressure becomes a flux-function. 

In addition, toroidicity affects qualitatively the eigenvalues and the transition points. 
We examined this impact by varying the aspect ratio a and found the following results: 

1. The eigenvalues Pq and F n ("compressible and incompressible) become lower as the 
aspect ratio decreases. For example, for a = 3 and a = 2 one finds Mq = 2.8 and 
Mq =2.3, respectively. 

2. The smaller a the lower the transition Mach numbers (Mo) m and A m (for any m). 
For example, for a = 3 and a = 2 the respective first-transition Mach numbers are 
(M )i = 1.692 and (M )i = 1.338 and the first-transition incompressible points are 
A x = -0.083 and A x = -0.448. 

The "compressible" transition velocities are found to be in general supersonic. Experi- 
mental supersonic toroidal velocities have been reported in Refs. ^1] and [13] . Owing 
to the above mentioned dependence of the transition points on a, however, it is possible 
to have subsonic transitions for appropriate low values of a. Thus, in this case the tran- 
sitions may be realized easier in spherical tokamaks; the minimum subsonic value of the 
first-transition point (M )i is 0.62 and corresponds to a compact toroid (a = 1). 

IV. Shafranov shift 

We evaluated the impact of the flow and the aspect ratio on the Shafranov shift of 
eigenstates with a single magnetic axis. (It is reminded that the Shafranov shift is defined 
as the displacement of the magnetic axis with respect to the geometrical center of the 
configuration, Ax = x m . a . — 1.) The results are summarized as follows. 

1. As Mq increases or A decreases the Shafranov shift increases. This increase for 
" compressible and incompressible equilibria can be seen in Tables Q and 121 respec- 
tively. It is interesting to note that for large positive values of A the Shafranov 
shift can become negative. As an example for a = 3 and A = 0.09 the shift is 
—0.0274. This means that the magnetic surfaces in this case are shifted inward. An 
inward shift of magnetic surfaces associated with poloidal flow in quasi-isodymanic 
equilibria was reported in Ref. [TSj. Also, suppression of the Shafranov shift by a 
properly shaped toroidal rotation profile was found in Ref. [T7j . 

2. The lower the aspect ratio the larger the Shafranov shift. This is shown in tables El 
and 0] for two different values of the pressure parameter P . 
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M 


Shafranov Shift 


0.1 


0.054 


0.4 


0.058 


0.6 


0.063 



Table 1: The Shafranov shift, Ax = x m , a . — 1, for various values of the Mach number 
(M ) and aspect ratio a = 3. 



A 


Shafranov Shift 


0.010 


0.045 


0.006 


0.049 


-0.001 


0.055 



Table 2: The Shafranov shift, Ax = s m a — 1, for various values of the incompressible-flow 
parameter A and aspect ratio a = 3. 



Aspect ratio 


Shafranov shift 


3 


0.092 


2 


0.150 


1.5 


0.209 



Table 3: The Shafarnov shift, Ax = x m , a , — 1, for P = 12 kPa and various values of the 
aspect ratio for the "compressible" case. 



Aspect ratio 


Shafranov shift 


3 


0.053 


2 


0.140 


1 


0.500 



Table 4: The Shafarnov shift, Ax = x m , a , — 1, for P = 110 kPa and various values of the 
aspect ratio for the incompressible case. 



V. Summary and Conclusions 

Equilibrium eigenstates of a magnetically confined plasma with toroidal flow sur- 
rounded by a boundary of rectangular cross-section have been investigated within the 
framework of ideal MHD theory. "Compressible" flows associated with uniform temper- 
ature but varying density on magnetic surfaces and incompressible ones with uniform 
density but varying temperature thereon have been examined on the basis of respective 
reduced equations [ (fTHj) and (|2*T|) ] and exact solutions [ (J23j) and (|23|) ]. The flow effect on 
the magnetic topology of the eigenstates has been examined by means of the parameters 
Mo and A associated with the quantities u 2 /T in the "compressible" and quj 2 in incom- 
pressible case, respectively. The exact "compressible" solutions considered are shearless 
[(uj 2 /(T))' = 0] while the incompressible ones have non-zero shear [(quj 2 )' ^ 0]. 

Owing to the flow one can consider either pressure eigenvalues, (Po) n n — 1)2,... 
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with the flow parameter F being free (F represents either Mq or A) or alternatively flow 
eigenvalues F n with free P . For fixed F in the former case and fixed P in the latter 
one, the eigenvalues satisfy the relations P n+1 > Pq and F n+1 > F n , respectively. The 
respective eigenfunctions for the poloidal magnetic flux-function tp can describe multi- 
toroidal configurations with n magnetic axes located on the mid-plane z — 0. When M 
is increased or A is decreased continuously there are transition points F m (m — 1, 2, . . .) 
at which an additional magnetic axis appears. Alternatively, by varying continuously the 
pressure parameter P there are transition points (Po) m associated with flow eigenvalues at 
which an additional magnetic axis is formed. This change in magnetic topology, possible 
only in the presence of flow, is crucially related to toroidicity because in the limit of in- 
finite aspect ratio the equilibrium equations are flow independent. The above mentioned 
"flow-triggered" transitions can be approximately viewed as a quasi-static "evolution" 
of the plasma by continuous flow variation through pressure eigenstates or alternatively 
by continuous pressure variation through flow eigenstates. The transition points have 
the following dependence on the aspect ratio a: the lower a the smaller (M ) m in the 
"compressible" case and A m in the incompressible one. 

Also, we have examined the impact of the flow and the aspect ratio a (as a measure 
of the toroidicity) on the Shafranov shift. The results show that the shift (a) increases 
as M take larger values and A smaller ones and (b) increases as a takes lower values. 
Furthermore, for large positive values of A the shift can become negative. 
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